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ABSTRACT 

We describe a family of circular, and elliptical, finite disks with a disk potential 
that is a power of the radius. These are all flattenned ellipsoids, obtained by squashing 
finite spheres with a power-law density distribution, cutoff at some radius R Q . First 
we discuss circular disks whose circular rotation speed v oc r°, with any a > —1/2. 
The surface-density of the disks is expressed in terms of hypergeometric functions 
of the radius: £ Q (r) oc (r/Ro) 2 ^ [1 - (r/R f} 1/2 ^(1 - a, 1/2; 3/2; 1 - R 2 Jr 2 ). 
The disk then has the same rotation curve as the sphere within the cutoff radius 
(up to a constant factor). We give closed expressions for the full 3-D potentials in 
terms of hypergeometric functions of two variables. We express the potential and 
acceleration in the plane at r > R Q , and along the rotation axis, in terms of simple 
hypergeometric functions. All the multipoles of the disk are given. We then generalize 
to non-axisymmetric disks with midplane axis ratio p < 1. The potential in the 
midplane is given by ip(r, 4>) oc r 2a F 1 (a + 1/2, —a, a + 1/2, a + 1; r sin 2 0; r) + c, where 
F ± is the hypergeometric function of two variables, and r = 1 — p 2 . For integer values 
of 2a the above quantities are given in more elementary terms. All these results follow 
straightforwardly from formulae we derive for the general, cutoff, power-law, triaxial 
ellipsoid. 



Subject headings: Disk galaxies, potentials 



1. Introduction 



Models of disk galaxies with explicit expressions for the disk potential and for the surface 
density are quite useful in studies of galactic structure, for example as initial states for studies 
of stability, spiral-arm formation etc.. Infinite-extent, power-law disk models-with power-law 
rotation curves: v(r) oc r a , and power-law surface densities: S(r) oc r 2a_1 — have their uses in 
this context (e.g Binney and Tremaine 1987 , Evans and de Zeeuw 1992 , Evans 1994 ). They are, 
however, rather unrealistic. Qian 1992 describes an even wider class of disks with known surface 
densities and potentials in the plane of the disk; these too are of infinite extent. We hope to 
contribute to the subject by describing a family of finite disks with power-law potentials (within 
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the material disk). These are produced by squashing spheres with a power-law density profile, 
cutoff at some radius R Q , into ellipsoids (disks in the limiting case), the rotation curve in the 
equatorial plane, within R , is then the same as that of the sphere (up to a constant factor). 
I 1 

Milgrom 1989 used this method to determine the surface density of a constant-acceleration 
(isodynamic) disk with v oc r 1/2 . We describe here the general case since, as far as we know, this 
useful class of models does not seem to be part of the known lore. In the Appendix we derive 
expressions for the gravitational field of a cutoff, power-law, triaxial ellipsoid. Prom these follow, 
after some straightforward manipulations, all the results of the paper. In section 2 we discuss thin, 
and thick axisymmetric disks; in section 3 we derive their full 3-D potentials; and in section 4 we 
describe the generalization to non-axisymmetric disks, which may also model galactic bars. 



2. Axisymmetric disk models 

Start with a spherical mass distribution, with density 

P = AR\ (1) 

cutoff at R = R Q . Now squash the configuration in one direction, by a constant factor q < 1, to 
obtain an ellipsoidal mass distribution. The gravitational acceleration, g(r), at a point in the 
midplane, a distance r < R Q from the centre, in the equatorial plane, is determined only by the 
mass interior to the ellipsoid going through r (by "Newton's third theorem"; see e.g. [Binney and 



III 

emaine 1987|). In particular, g(r) is independent of the cutoff radius R Q (provided R Q > r). The 



only dimensional parameter characterizing the system, on which g(r) may depend is then A, to 
which, in fact, g must be proportional. On dimensional grounds, the only form that g(r) can have 
is then 

g(r) = -QAGr J+1 , (2) 

where Q is a dimensionless constant that depends on the axis ratio and on 7. This argument 
generalizes to non-axisymmetric disks, and indeed to the interior field of a cutoff, power-law, 
triaxial ellipsoid (see section 4 and the Appendix). We see that the equatorial rotation curves 
of this family of ellipsoidal mass distributions is a power-law within the cutoff radius, with 
v(r) = (QAG) 1/2 r a , a = 1 + 7/2. Beyond R Q , g(r) does depend on R a , and the above dimensional 
argument breaks down, as the expression for g(r) may now contain a factor f(r/R Q ). 



The surface density is just that of the spherical distribution 

Y, a {r) = 2A (r 2 + x 2 y /2 dx. (3) 



Defining u = r/R Q , s = [(R /r) 2 - 1] 1/2 = (u~ 2 - 1) 1/2 , and t = 2ARl+\ we can write after a 
change of variable 

£» = is^ +1 s /' A- 1/2 (l + s 2 A)^ /2 dX. (4) 
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This, in turn, can be written in terms of the hyper geometric function 

£<*(r) = tu 2(a ^(l - u 2 ) 1 ' 2 2 F X (1 - a, 1/2; 3/2; -s 2 ), (5) 



(see Gradshteyn and Ryzhik 198C 3.197.3). We see that near the disk's edge, where s ~ 0, and 
2 F 1 « 1, S Q (r) « S(l - u 2 ) 1 ' 2 , for all the disks. 

For integer values of 2a the integrals in eqs.(||)(||) can be expressed in terms of elementary 
functions. Some interesting special cases are: a = 1-which in the thin-disk limit gives the Kalnajs 
disk ( Kalnajs 1972| )-has 

E(r) = S(l--u 2 ) 1/2 ; (6) 
a = 0-which is the Mestel disk ( Mestel 1963| ) in the same limit-has 

£(r) = tu- 1 tan- 1 ^" 2 - 1) 1/2 ; (7) 

and the isodynamic disk ( [Milgrom 1989 ) with constant acceleration (a = 1/2) has 

£(r) = Ssinh- 1 ^- 2 - 1) 1/2 . (8) 

The limit a — > —1/2 corresponds to a point mass at the origin (a Keplerian rotation curve), and 
smaller values of a are excluded as they correspond to an infinite total mass. 

For a a positive integer, 2 Fi m eq.@ is a polynomial of order a — 1 in s 2 , and so, S(r) is, up 
to the factor (1 — u 2 ) 1/2 , a polynomial of (even) order 2(a — 1) in r. 

Values of a > 1 correspond to disks with non-monotonic surface densities: for 1 < a < 3/2 
S(r) — S(0) oc r 2(t * -1) near the origin; for a = 3/2 S(r) — S(0) cx — r 2 lnr there; and for a > 3/2 
S(r) — S(0) cx r 2 . For all these cases £ is then increasing near the origin, but decreases further 
out. 

The proportionality constant, Q, in expression (|2|) for g may be obtained by specializing 
results from Appendix A. 

Q(k, a) = — — ^(1/2, a + 1/2; a + 3/2; «), (9) 
2a + 1 

where k = 1 — q 2 . 

The interior potential for < i? [l + (n/q 2 ) cos 2 9]~ 1/2 (for 7 7^ —2) as a function of the polar 
coordinates R and 9 is 

wtH, 9) = 27TAG , R 2a+1 F 1 (a + 1/2, -a, a + 1/2, a + 3/2; /tsin 2 fl, k) + c, (10) 
a(2a + 1) 

Where is the first Appell function, or the first hyper geometric function of two variables (see 
Gradshteyn and Ryzhik T980| 3.211), and c = -27ra" 1 AGR 2 Q a k~ 1/2 sin" 1 k 1 ' 2 . 
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Hereafter we concentrate on infinitely thin disks, for which the results may be put in simpler 
terms. For these 

Q = Q {a) = 27r 3/2 r(a + l/2)/r(a + 1). (11) 
The potential in the disk is given (except for the Mestel case) by 

^ (r) = ^- + C - (12) 
The constant c can be determined from the potential outside the disk treated in section 3: 

ir 2 AGRl a MGf(a + l/2) , , 

c= = — p 5 , (13) 

where 

M = AnAR 2a+1) /(2a + 1) (14) 
is the total mass of the disk. Note that <p, unlike g, does depend on R , in general. 
The total energy on a circular orbit at radius r < R Q may now be written 

r \ 2a (a + l)r(a + l/2)" 

T ) r(i/2)r( a + i) . 

This energy is negative everywhere in the disk for a < 1, vanishes just at the edge, for the 
Kalnajs case, and is positive in some region near the edge for a > 1, viz. for r > rj(a)R , where 
77(a) = [( a + l)r(a + l/2)/7r 1/2 r(a + 1)]~ 1/2q . It can be shown that 77 < 1 for all a > 1, and we 
find numerically that it has a minimum of ~ 0.96, at a ~ 4.37; it is thus only in a narrow region 
near the edge that E(r) > 0. We have seen that a > 1 disks have nonmonotonic surface densities. 
Such disks may not be so realistic models by themselves; they are useful as components when we 
expand the potential of more general disks. We may find advantage, e.g., in disks whose potential 
is a finite sum of power laws, some possibly with a > 1. 

As an example, consider two of our disks with a 1 < a 2 , but with the same R Q , and the same 
S = IARq"' 1 (the values of A are different). The corresponding, unsquashed spheres then have 
the same density at R Q . We see from expression (|3|) for the surface density that, for a fixed S, 
S a (r) is a monotonically decreasing function of a for all r < R Q . Thus 

E" 1 " 2 = E" 1 - E" 2 (16) 

is positive everywhere. The acceleration in such a disk is 



E(r) = MG <a + V2) 
{ ' 2R D a 



1 



(15) 



g = -tGir [i/2 



r(q 1 + l/2) ^ 2ai _ 1 _ ^ + 1/2) ^^ 



(17) 



r(«! + 1) r(a 2 + 1 

(u = r/R ), and is directed inward everywhere, because u 2a ~ 1 T(a + l/2)/T(a + 1) is decreasing 
with a. Furthermore, S Q1Q2 (r) vanishes at the edge faster than (1 — « 2 ) 1/2 ; in fact, there 

S^(r)«^E(a 2 -aO(l-n 2 ) 3/2 , (18) 
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which can be seen by taking 2 Fi to first order in s 2 in eq.(|5|). We can get disks that vanish even 
faster at the edge by taking the difference of two such disks with different pairs of a's with the 
same difference. In the limit a 2 — > ot-i we get a disk whose surface density is the derivative of E a (r) 
with respect to a at a fixed S. 



3. The field outside the axisymmetric thin disk 



Building on formulae in Appendix A we get for the potential at a general position outside the 
disk (for 7 7^ — 2-the Mestel disk requires separate treatment) 



cp(r, z) 



MG 
R„ 



(p(r/R ,z/R ), 



(19) 



where the dimensionless potential <p, as a function of the dimensionless variables 
g = r/R , ( = z/Ro, is 



<p(q, 



a + 1/2 



a 



sin-^l + L) 



-1/2 



1 r w t«-va(A a - (H) a 



2 J 



[1 - ty- 1 



/2 



(It 



Here, A is the dimensionless polar radius A 2 = £ 2 + g 2 , 

L = ^{A 2 - 1 + [(A 2 - l) 2 + 4C 2 ] 1/2 }, 



(20) 



(21) 



(L is lj(1,R) of appendix A, with Y = R a ), and W = 1/(1 + L). The integral in eq.(20) may be 
identified, after some manipulations, as an integral representation of the first Appell function. We 
thus get finally 



<p(q, C) 



a + 1/2 



'1/2 



a 



where 



A 2Q W^ +1/2 i?(a + 1/2, - a , a + 1/2, a + 3/2; U, W), 

2a 



U = 1 - L/A 2 . 



(22) 
(23) 



For all non-zero integer values of 2a the integral in eq.(|20|) can be performed and expressed in 
terms of elementary functions. One case in point is the Kalnajs disk which is just an homogeneous 
ellipsoid, for which the potential is given in McMillan 1958| or in Binney and Tremaine 1987 . 
Another instance is that of the constant-acceleration disk where, integrating explicitly, we get for 
the dimensionless 3-D potential 



<p(g, C) = 2L 1/2 - 2A + 2 edn-^l + L)- 1 ' 2 + |C| In 



(LV 2 -|C|)(A + |C1) 
(LV 2 + | C |)(A-|C|) 



(24) 
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For the Mestel case with a = 0, eqs.(20) and (22) give 0/0, and it requires special treatment, which 



does not seem to yield a simple closed expression. In general, all expressions for the accelerations 
hold for the a = case; those for potentials and energies require either a separate derivation, or 
the application of the L'Hospital rule. For example, we can use this rule with eq.(|2"2|) to obtain for 
the Mestel case 

<p(Q,0 = -^(lnA^lnWOsin-^l + L)- 1 / 2 - 
1 d 

~2 Wl ' 2 d^ Fl{a + 1/2, ~ a ' a + 1/2, a + 3/2; u ' w) ^=°- (25) 

The integral in eq. (p0|) , for can also be expressed in terms of simpler functions when 

we confine ourselves to the plane (z = 0, r > R ), or to the symmetry axis (r = 0). We then get 
for the dimensionless potential for the former case 

${Q > 1, 0) = ^±1 sin" 1 q- 1 - J- 2^(1/2, a + 1/2; a + 3/2; ^ 2 ); (26) 
2a lag 

along the axis we have 

m C) = 2 A^- sin-^l + CY 1 ' 2 - ^(1, « + 1/2; « + 3/2; -C 2 ). (27) 



[Equations (Pq) (|27|) can also be obtained as special cases of eq. ([22]) ; for £ = we have U = W, and 
for g = 0we have U = 0; in both cases F 1 can be expressed as a hyper geometric function of one 
variable.] 

The acceleration in the plane can be brought to the form (valid also for a = 0) 

g(r >R o ,z = 0) = a.FUl/2, a + 1/2; a + 3/2; R 2 Jr 2 ). (28) 

For r — > oo 2 -Fi - * 1, and we get the Keplerian expression. At the disk's edge 
2^(1/2,0 + 1/2; a + 3/2; 1) = vr 1/2 (a + l/2)r(a + l/2)/r(a + 1) is the ratio of the disk's 
field to that of a sphere with the same mass 

The acceleration along the symmetry axis can be written as 

g (r = 0,z) = -^- 2 F 1 (l,a + 1/2; a + 3/2; -R 2 Jz 2 ). (29) 



Using relations found in Magnus and Oberhettinger 1949 to get explicit expressions for 2 F 1 



in special cases we obtain from eqs.(^) and ( p9[ ) (with t = g 1 = R /r, and t = C 1 = Ro/z): 
for the Kalnajs disk 

g{r > R0 , z = a) ^ ^^-'^-^ (30) 
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(reproducing a result of McMillan 1958| ), and 



MG 3 / tan" 1 A . , 

9(r = , ,)=-_-(! — ). (31) 

For the isodynamic disk 

g(r > R , z = 0) = ~ — , (32) 

and 

MG ln(l + t 2 ) 

5 (r = 0,z) = . (33) 

z 1 t 1 



For the Mestel disk we get the Mestel 1963| result for the plane 

MG « , , <. 

g(r > R o ,z = 0) = -—t- 1 sm- 1 t, (34) 

and 

MG tan" 1 1 , s 

g(r = 0,z) = -— — , (35) 

along the axis. 

Because 2 F\ is a power series with simple expressions for the coefficients, any of 



eqs,(p6|)(27)([2q) (p9[) gives straightforwardly the multipoles of the disk: Write the potential (using 
polar coordinates) as 

oo 

ip(R, 9,<I>) = J2 n e R- (e+1) P e (cos 6) (36) 

1=0 

(only even £'a appear). Use, e.g., eq.(|2~8|); the hypergeometric function appearing in it is a power 
series 

oo 

^(1/2, a + 1/2; a + 3/2; R 2 Jr 2 ) = ]T f k R 2 k r~ 2k , (37) 

with 

= r(fc + l/2)(q + l/2) 
Jk TT^ 2 k\(a + k + l/2) " 1 ' 

The radius of convergence of this series in R /r is 1 (i.e. the expansion is valid for r > R ). 
Comparing with the R derivative of expression (36), putting P 2 fe(0) = (— l) fc (2/c)!2~ 2fc (A;!)~ 2 , and 
using some identities for the T functions, we get for the multipole coefficients of the disk 

^ = (-ifGMRf — " + ^ — . (39) 

We see that, up to a numerical factor, jjL 2k is given by GMR 2 ^, which is to be expected on dimensional 
grounds as, in constructing fx, we can only avail ourselves of the dimensional quantities M and R D 
(or A and R Q ). For example, the quadrupole moment is fj, 2 = —GMR 2 Q (a + l/2)/3(a + 3/2). We 
reiterate that while expressions (]2"0D fl22]) for the potential are valid everywhere outside the disk, 
the multipole expansion is valid only outside the sphere of radius R Q . 
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4. Non-axisymmetric thin disk models 



Our disks may be generalized to non-axisymmetric ones, which may serve as models for 
non-axisymmetric galactic disks or for bars. Start again with our cut-off-power-law sphere and 
squash it by a factor p along the x axis and into a thin disk along the z axis, to obtain a finite, 
elliptical disk of axis ratio p. The field vanishes within a triaxial homoeoid obtained by squashing 
a spherical shell of constant surface density by different factors along different axes. Thus, our 
dimensional argument of section 1 still holds, and the acceleration within the disk itself must still 
be proportional to r 7+1 , and oblivious to the cutoff. Thus 

g r (r, 4>) = -AGr-< +1 P r (T,"f, sin 2 0), (40) 

g+(r, <j>) = -AGr^P^r, 7, sin 2 0), (41) 

where r, <j> are polar coordinates, and r = 1 — p 2 is a measure of the departure from axisymmetry. 
When 7 ^ —2 eq.(^0|) tells us that the midplane potential inside the disk must be of the form 

AGr 1+2 

V(r, <P) = -—^P r (r, 7 , sm 2 <P) + C(</>, R Q ) (42) 
7 + 2 

(note that the potential, unlike the acceleration, does depend on R ). Deriving g^ = — r_1 §^ from 
eq.(f42|) and comparing with eq.(fi"I|) we see that C may not depend on <j) and that 

All that remains is to derive one angular factor P r . For the special case 7 = —2 eq.(^) tells us 
that the potential must be of the form 

tp(r,<f>) = AG In rP r (r, 7, sin 2 <j>) + C(<f>, R a ), (44) 

but now for eq.@ to hold we must have ^ = 0, C = AGC(4>) + c(R Q ), and P* = We learn 
then that for this case the radial acceleration is (ft- independent. 

The surface-density of the non-axisymmetric disk is obtained simply from expression (||) for 
the axisymmetric case by multiplying the latter by a factor 1/p-to account for the squashing 
along x-and by replacing the variable r in eq.(|5|) by r[l + (r/p 2 ) cos 2 0] 1/2 ; the boundary is at 
r = R /[l + (r/p 2 )cos 2 0] 1/2 . 

We have derived the field inside the disk by integrating the contributions of thin triaxial 
homoeoids to the exterior potential (taken from table 2-1 of Binney and Tremaine 1987) weighted 
by the power-law density profile. The derivation is given in Appendix A. The components of the 
acceleration may be written as 

,~ a 1 r t a - 1/2 (l - Tsin 2 (/> t) a , 
g r = -2TiAGr 2a - x / )—, dt, (45) 
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g. = 27rAGr 2Q - 1 rsin0 cos0 / , , / \ 1 , dt, (46) 

where a = 1 + 7/2. The integrals in eqs,(45)(|46|) may be written in terms of the first Appell 
function 

g r = -QoAGr^F^a + 1/2, -a, a + 1/2, a + l;r sin 2 <£; r), (47) 

ff0 = Q AGr 2 "- 1 7 + , r sin cos x 
(a + 1) 

xF 1 (a + 3/2, -a + 1, a + 1/2, a + 2; r sin 2 0; r). (48) 

These are valid for all values of 7 > —3, including 7 = —2 (a = 0). In the last case we see from 
eq.((45|) that indeed g r does not depend on </>. 



When 7 7^ —2, the potential in the disk is given by 
O AGr 2a 

<p(r,4>) = - £2 -= Fx (a + 1/2, -q, a + 1/2, a + 1; r sin 2 0;r) + c, (49) 

2a 

where 

c= 2_ 3jFi l/2,l/2;l;r 50) 

a 

is calculated in Appendix A. In the axisymmetric case 2 Fi = F 1 = 1 as the arguments vanish and 
we get the results of section 2. 

For 7 = —2 the interior potential is 

<p(r,<l>) = 2n 2 AG 2 F 1 (l/2,l/2;l;r)ln(r/R ) + 
[i ln[(l -r S m^t)/t] 

-vAgC l H{l ; TtW) \ l dt, (51) 
Jo t 1 / 2 (l -t) 1 /2(l -rt) 1 / 2 ' v ; 

to be compared with eq.(44). The first term gives g r , the second g^, and the third contributes to 
the constant that corresponds to ip = at infinity. 

We next discuss various special cases. As in the axisymmetric case non-zero integer values of 
2a afford simplifications. We get for a = 1: 

V?(r,0) = K 2 AGr 2 [ a F x (3/2, 3/2; 2; r) - ^(3/2, 5/2; 3; r)] + c, (52) 

which is simply an anisotropic, harmonic potential (known to be the case inside an homogeneous 
ellipsoid-it is the basis, e.g., of the Freeman bar model described by Binney and Tremaine 1987 ). 

For the a = 1/2 disk we have 

<p(r, (f>) = 4nAGrB(T, sin 2 <f>) + c, (53) 
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where c is given by eq. (|50|) , and 

„ I cos 01 , / t cos 2 d>\ 1 I sin 01 1 + r 1/2 | sin 01 . 

5 = - — ! — —P — -tan -1 +- P-ln -f. (54) 

(l_ r )i/2 T i/2 \ i_ T J 2^/2 1- raisin 0| v ; 

The weak-asymmetry limit, r <C 1, is readily obtained for the a general a, as F x is a known 
power expansion in its two arguments. To first order in r we have 

Q AGr 2a r , r(a + l/2) ^ , o ^ h r ,,,, 
y(r,0) = — ^— [1+ 2(a+1) (l + 2acos 0)] — [I+4L (55) 

5r = -QoAGr 2 - 1 !! + ^fe^r^a + 2a cos 2 0)], (56) 

2(a + 1) 

and 

<7 = Q ^Gr 2 - 1 ^^ r sin cos 0. (57) 

The disk potential and (radial) acceleration along the principal axes may also be written in 
simpler terms for general a. Along the (shorter) x-axis, where = 

g r = -Q AGr 2a - x 2 F Y {a + 1/2, a + 1/2; a + 1; r), (58) 

and along the longer axis 

g r = -QaAGr 2 "- 1 2 F 1 (l/2, a + 1/2; a + 1; r). (59) 

The most general case pertains to the external field of a thick disk. This is given in 
Appendix A in terms of a Lauricella function of four variables, F^\ which is a generalization 
of the hypergeometric function (see Exton 1976] ), In any of the three symmetry planes of the 



ellipsoid this reduces to a Lauricella function of three variables F^ (F^ is the same as F^. It is 
straightforward to obtain all these from formulae in Appendix A. All these functions are single 



integrals between and 1 of simple algebraic functions. [Exton 1976j gives a short FORTRAN 
program to evaluate these functions. 



We thank David Earn for useful comments. 



A. Formulae for the field of a cutoff, power-law ellipsoid 



We start with the expression, given in Table 2-1 of Binney and Tremaine 1987 for the potential 



at position R outside a thin, triaxial homoeoid of axes a i; 1 = 1, 2, 3 

G 



<p(R) 



-M shell K{b), 



(Al) 
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where 6; are the axes of the ellipsoid that is confocal with the shell and passes through R: 
&? = »? + A (a, R), and A is defined by 

£l*/(oJ + A) = l. (A2) 

The quantity K is given by 

K{b)= 'o m+sM+sM+s)]^- (A3) 

We now calculate the potential outside a cutoff, power-law ellipsoid, whose external axes are pY, 
Y, and qY, integrating the contributions of thin-shells of axes kpY, kY, kqY, with < k < 1. The 
mass of the shell between k and k + dk is that of the original spherical shell dM = 4nAY'' +:i k' 1+2 dk, 
so that 

<p(R) = -2ttAGY-< +2 J, (A4) 

with 

j — Y b 1+2 dk I — ( A^) 

J J {[(pkYy + x + s ][(kYY + x + s ][(qkYy + x + s]y/ 2 ' { ' 

Changing variables to t = (A + s)/(&y) 2 we have 

rl roo 

J = / W +1 dk / /(t) dt, (A6) 

with 

/(i) = [(/ + i)(l + i)(<Z 2 + i)]- 1/2 (A7) 

that is independent of R, k, or Y, the dependence on which enters through the lower integration 
boundary lo = X(k, R)/(kY) 2 . We now integrate by parts over k (the case 7 = —2 requires a 
special treatment) to get 

1 f°° 1 Z* 1 di ) 

J = — - \ f(t) dt + — dk k^f(u). (A8) 
7 + 2 Ju,(\,A) 7 + 2 Jo dk 

In the second integral-call it J- we change variables to u, noting that, from the definition of A in 
eq. ( |A2|) we have 

f 2 v 2 C 2 

p 2 - + u 1 + u q + to 

Here £ = x/Y, v = y/Y, ( = z/Y. At k = u(0, R) = 00. We then get 

r°° f e v 2 c x 2 ^ 

J = - dco f(u J )(^— + -— + ^—) . (A10) 

Juj(1,R) \p 2 +UJ 1+lu q 2 +ujj 

It is useful to change variables to t = 1/(1 + u), and defining r = 1 — p 2 , k = 1 — q 2 we can write 
j = f w t h+1)/2 [g 2 (l - Kt) + ^ 2 (1 - Kt)(l - rt) + C 2 (l - rt)p+ 2 )/ 2 

X [(1-Kt)(l-Tt)](^+ 3 )/ 2 ' 1 j 
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where W = 1/[1 + uj(1,R)]. For R on the surface of the cutoff ellipsoid lo(1,R) vanishes, so the 
upper integration boundary is 1. 

To obtain the acceleration at position R inside the ellipsoid we make use of the fact that the 
layers outside R do not contribute. Thus we can use the above expression for the potential outside 
an ellipsoid that is cutoff beyond that through R, then take the R gradient, and calculate the 
integrals at the boundary. When taking the R gradient, we note that the two contributions from 
differentiating the boundaries of the integrals cancel, and we are left with the following expression 
for the acceleration inside the cutoff ellipsoid 

g = -^^V[R-' +2 I(x,y,z)], (A12) 
7 + 2 

with I a function of the polar angles, but not of R. 

_ t^ 2 [ X 2 (l - Kt) + £ 2 (1 ~ Kt)(l - Tt) + Z 2 (l - Tt)p+ 2 >/ 2 

1 ~ Jo [(l-«t)(l-rt)]^ a ' ( } 

where x = sin#cos</>, y = sin 6 sin 0, and z = cos 6 are the angular factors of the coordinates R. 
This rather restricted form of the field follows from the same dimensional arguments we applied 
in sections 1 and 4. The interior potential is then (7 ^ —2) 

p _^r s /(i,P)H (A14) 
7 + 2 

By comparison with the exterior potential at the boundary, where W = 1 we get 
2irAGRZ+ 2 s , 4ttAGR 



r°° 4ttAGR 1+2 

/ f(t) dt = - ^ U ° fi(l/2, 1/2, 1/2, 3/2; k,t). (A15) 
7 + 2 Jo 7 + 2 

To get our results in the midplane of a disk we take 9 = tt/2, and go to the limit q = (k = 1). 
All our expressions for the interior acceleration then follow. Our expression for the potential outside 
an axisymmetric disk eqs.([T^) (pC|) is obtained by putting r = and Y = R Q in eqs. (|A4|) (|A8| ) ( All ), 
and then writing the resulting integrals in terms of hypergeometric functions. 



Regarding the field at an arbitrary position outside the ellipsoid, we note that in eq.( All| ) 



the term in square brackets in the numerator of the integrand may be written in the form 
(R/R ) 2 (l — bt)(l — ct). Thus, by changing variables to t = sW we can write 

j = -B f s wl)/2 {{\ - bWs)(l - cVFs)p +2)/2 [(l - kWs){1 - tWs)Y { ~i +3)/2 ds, (A16) 
Jo 

with B = [R/ R ) 2a W a+1/2 . This can be written in terms of the Lauricella function of type D, 
which is a generalization of the hypergeometric function to many variables (see |Exton 1976 2.3.6). 



The first integral in eq.(A8) can also be brought to a simpler form, and we finally get for the 
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potential outside the cutoff, power-law ellipsoid 



cp(R) 




~W a+1/ \a + \/2)- 1 x 



-{2W 1/2 F 1 (l/2, 1/2, 1/2, 3/2; kW,tW) - 



xFg>(a + 1/2, -a, -a, a + 1/2, a + 1/2; a + 3/2; bW, cW, kW, tW)}. (A17) 



We remind the reader that W = (1 + co) 1 , where w is the positive solution of 



The exterior potential may be written in terms of hyper geometric functions of three variables 
when we are in any of the symmetry planes of the ellipsoid. 
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